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ABSTRACT 

The temporal-spectral evolution of the prompt emission of gamma-ray bursts 
(GRBs) is simulated numerically for both leptonic and hadronic models. For weak 
enough magnetic fields, leptonic models can reproduce the few seconds delay 
of the onset of GeV photon emission observed by Fermi-LAT , due to the slow 
growth of the target photon field for inverse Compton scattering. For stronger 
magnetic fields, the GeV delay can be explained with hadronic models, due to the 
long acceleration timescale of protons and the continuous photopion production 
after the end of the particle injection. While the FWHMs of the MeV and 
GeV lightcurves are almost the same in one- zone leptonic models, the FWHM of 
the 1-30 GeV lightcurves in hadronic models are significantly wider than those 
of the 0.1-1 MeV lightcurves. The amount of the GeV delay depends on the 
importance of the Klein-Nishina effect in both the leptonic and hadronic models. 
In our examples of hadronic models the energies of the escaped neutrons are 
comparable to the gamma-ray energy, although their contribution to the ultra 
high-energy cosmic rays is still subdominant. The resulting neutrino spectra 
are hard enough to avoid the flux limit constraint from IceCube. The delay of 
the neutrino emission onset is up to several times longer than the corresponding 
delay of the GeV photon emission onset. The quantitative differences in the 
lightcurves for various models may be further tested with future atmospheric 
Cherenkov telescopes whose effective area is larger than that of Fermi-LAT , such 
as CTA. 
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radiation mechanisms: non-thermal 
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Introduction 



The nature of the prompt emission mechanism of gamma-ray bursts (GRBs) is still con- 
troversial, and remains a challenging problem in high-energy astrophysics. Most of the GRB 
spect ra peak around th e 0.1-1 MeV range, and on the whole are well fitted by the Band func- 



models have been proposed such as photospheric em ission (see e.g. iRees fc Meszarosl 12005 



' Piran 


2005; 


Meszaros 


3. However, alternative 


Rees & Meszaros 


2005; 



Giannios 1120071 ; lBeloborodovll2010t iRyde et al. II2010I ). The recent detection of GeV photons 
with Fermi-LAT has opened up the possibility of constraining such models. In some objects 



Fermi has also found in the GeV energy range addit i onal spectral component s flAbdo et al 



2009a; Ackermann et al. 2010; Abdo et al. 2009b; Ackermann et al 



201 1 ). Moreover, a 



general feature is that the onset of the GeV emission tends to be delayed relative to the 



onset of the main MeV emission (I Abdo et al. I l2009d . and others including the above ref- 
erences). The typical delay timescales for the lon g Fermi-LAT G RBs are 3-5 s as seen in 



GRB 080916C, GRB 090902B and GRB 090926A flMeszaros 1120121 ) . The simplest model to 
explain th is would be i nverse Compton emission from t he same internal dissipation region as 
the MeV ( Corsi et al. 2010 ; Asano fc Meszaros 12011, hereafter AM 11). A hadronic origin 



for th e GeV emission have also been proposed (lAsano et al. Il2009al . |2010| ; iRazzaque et al. 
20101 ) . although such models require a much larger energy in accelerated protons than in 
the emitted gamma-ray energy. Alternatively, if the Lorentz factor of the external shock 
propagating in the interstella r medium is as high as 1000, the delayed onset of th e GeV 
emissions may be attributed (jGhisellini et al. 1 12010| ; Kumar fc Barniol Duran 1 12010| ) to an 
early onset of a forward shock synchrotron afterglo w. Other mechanisms to repro duce the 
GeV delay are proposed by several authors (see e.g. iToma et al. Il2009i ; lloka 1120101 ). but we 
do not arrive at an established interpretation yet. 

The photon statistics above the GeV range provided by Fermi are not sufficient to 
distinguish between the internal or external origin of the high-energy emission. However, 
it is expected that future multi-GeV ob s ervations with atmo s pheric Cherenkoy telescope 
arrays such as CTA (jGilmore et al. 2012 ; Kakuwa et al. 2012 ; Inoue et al. 2012 ) will dras- 



tically improve the data quality, owing to their large effective area. Lightcurves and spec- 
tral evolution determinations expected from such telescopes should provide critical informa- 
tion on the GRB physics, currently being pioneered by Fermi. To discriminate between 



(e.g. 


Pe'er & Waxman 


2005 




Pe'er 


2008 


Bosnjak et al. 


2009; 


Daigne et al. 


2011) 



Belmont et al. 2008; Vurm fc Poutanen 2009 
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In this paper, we concentrate on the internal dissipation regions as possible models for 
explaining the delayed onset of the GeV emission. In AM11, we developed a time-dependent 
code that can follow the evolution of the particle energy distributions in a relativistically 
expanding shell with Lorentz factor T and initial radius Rq from the central engine. In this 
previous AM 11 calculation we included only leptonic processes. However, GRBs have been 



considered a s possible sources of ultra high-energy cosmic rays (UHECRs, IWaxman I Il995 



Vietri Ill995l ) , and the electromagn etic cascades initiated by photopion production from accel- 
erated protons (hadronic cascade, iBottcher &: Dermer Ill998l ; iGupta fc Zhang 1 120071 ) result 
in both neutrinos and very high energy photons. The latter, potentially, could reproduce the 
delayed onset of the GeV emission observed, because of the long proton-acceleration timescale 
and/or the continuous pion production even after the end of particle inj ection. Also, the flat 
spectrum of the keV-GeV power-law compon ent seen in GRB 090 902B ( lAbdo et al. Il2009bf ) 
is compatible with a hadronic cascade model fjAsano et al. Il2010f ). Thus, in order to explore 
the time evolution of hadronic models, in this paper we incorporate the full details of the 
hadronic processes into the time- dependent code from AM11. With this expanded code, we 
are able to calculate the time evolution of both the usual leptonic model, as well as the time 
evolution of the radiation from hadronic cascade models. 

The observed lightcurves can show several pulses, and the origin of the time variability 
timescale is unknown. One possibility for this is variability in the central engine ejection, in 
which case the observed pulse widths are limited by the usual formula ~ Ro/cT 2 , and there 
could be multiple MeV-photon sources before the GeV onset. In this case, the GeV-delay is 
due to a superposition of the effects from various shell sources with different model parame- 
ters, the first pulse tending to have a soft spectrum in this scenario. Alternatively, the vari- 
ability migh t be due other effects, such as hydrodynamics ! turbulence in the sh e ll or behind 



the s hocks (IKumar &: Narayan I l2009t IZhang et al.l l2009t iMizuno et al. 1 1201 It llnoue et al. 



201 ll ) . with various unconstrained parameters. To avoid such model dependent uncertain- 
ties, in this paper we restrict ourselves to the study of a single primary pulse produced in 
a one-zone model, our code producing a smooth lightcurve, which can be considered as an 
envelope which averages a possible underlying variability with some smoothing. What we 
test is whether a single pulse in such a one-zone models can reproduce the observed GeV 
delays or not, with either leptonic or hadronic radiation mechanisms. The differences in the 
lightcurves for the various models explored in this paper provide clues which would be use- 
ful for discriminating between them, using future observations from atmospheric Cherenkov 
telescopes. 

In §2J we go into some of the details of the new numerical code for calculating the 
spectral evolution. The basic model features and our numerical results are presented in SJ3] 
for the leptonic and in §H for the hadronic models, respectively. A discussion and summary 
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of these results is given §5j 



Numerical Methods 



The numerical code in this paper is developed from the one- zone code in AM11. With 
this code we can simulate the temporal evolution of photon emission. The photon source in 
this code is a shell expanding with the Lorentz factor V from an initial radius R = Rq. This 
shell is responsible for one pulse, as part of a lightcurve which may have multiple pulses. The 
calculation of the photon production (for details see AM11) is carried out in the shell frame 
(hereafter, the quantities in this frame are denoted with primed characters). In this paper 
we express particle kinetic energies of electrons, protons, and photons as e e = (% — l)m e c 2 , 
£p = (7p ~ l)m p c 2 and s, respectively. 

The code used in AM11 can simulate the injection and cooling for electrons/positrons, 
and production, absorption and escape for photons with physical processes of (1) syn- 
chrotron, (2) Thomson or inverse Compton (IC) scattering (including the Klein-Nishina 
regime), (3) synchrotron self-absorption (SSA), (4) 77 pair production, (5) adiabatic cooling. 
In the present paper, a new feature is the incorporation of the hadroni c processes (based on 



the si m ulation methods i n the s eries of GRB s tudies of lAsano I (120051 ); lAsano fc Nagataki 



fl2006f ): lAsano fc Inoue I (120071 ): lAsano et al. I (|2009bJ)) into the code of AM11. For GRB 
prompt emission, the most important process in hadronic cascade is photopion production. 
The timescale of photopion production in the one- zone approximation is written as 

1 



t 



1-1 

P7 



d/i'(l - n')n'Jef 



where n^(e) is the energy distribution of the photon density, \x is the cosine of the photon 
incident angle, and K P1 is the inelasticity. We adopt experimental r esults for the cross 



sectio ns <r P7 for p-f — > nn + , pit , n-n- + ir°, and p-K +r K for e" < 2 GeV (see lAsano fc Nagataki 



2006 



for details), where e" is the photon energy in the proton rest frame. Multi-pion 
production due to high-energ y gamma-rays far above the A-resonance (e" ~ 300 MeV) is not 
important in our simulations (IMurase &: Nagataki II2006I ). because the standard synchrotron 
model generates a sufficiently soft photon spectrum with the typical index of —1.5. For 
the pion production by nj, we adopt the same cross sections as pj. The inelasticity is 
approximated by a conventional method as 



ml 



ml 



(2) 



where s is the invariant square of the total four-momentum of the p7 (717) system, and 

For the double-pion production, we approximate the inelasticity 



m p is the proton mass 
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by replacing wi t h 2m^ . Adopting the same Monte Carlo method as lAsano I (120051 ) or 
Asano fc Nagataki I ( 120061 ). we estimate energy loss of protons/neutrons and inject pions 
every time step following equation ([I]). 

For reference, we rewrite equation ([1]) for a simple case, in which a source with the bulk 
Lorentz factor V emits photons at a radius R with a power- law spectrum n'^e') oc e' a for 
e' < £p Cak and luminosity Ll (integrated below £ pe ak)- Here, e pea k corresponds to the spectral 
peak energy of the Band function. The result becomes 



l-a 



in 



(2 + a)L L 2 



,/-(cx+2) /_( Q+ l) 



'peak 



for 



,th 



7* 



300MeV 



'peak 



(3) 



(4) 



If we adopt the observationally typical index a = — 1, the timescale does not depend on the 
energy of proton. As is well known, the photomeson pr oduction efficiency f P7 = t 'It' is 



close to unity for fiducial param eter sets (t' exp = R/cT, IWaxman fc Bahcall 1119971 ; lAsano 
20051 ; iMurase fc Nagataki II2006I ). The theoretically expected value a = —1.5 in the standard 
synchrotron model leads to more efficient pion production for higher energy protons as 
/ P7 oc 7p 5 . Additionally, the hadronic cascade may produce an even softer spectrum. Our 
time-dependent simulation is a powerful tool to follow the pion production rate as the photon 
spectrum evolves. 

In terms o f the Thomson optical depth r T , the pp reaction efficiency is / pp ~ 0.05tt 
(IMurase 1 120081 ). Thus, we can neglect the pp-collisions in the optically thin cases in this 
paper. However, since our code is planned to apply also to optically thick cases in the nea r 
future, we have mounted the pp reaction on our code, as was done in IMurase et al. I (120121 ). 



Calculated with the numerical simulation kit Geant4 ( Agostinelli et al~l 20031 ). we provide 



tables for the cross section, inelasticity and pion-multiplicity for pp-collision. The procedure 
to follow proton cooling and pion injection is the same as that for jry-collision. 

The Bethe-Heitler (BH) pair production process (jry — > pe + e~) is also taken into ac- 
count. The cooli ng and spectral pair injectio n rates are calculated with the cross section and 



inelasticity from IChodorowski et al. I (119921 ). 



In our code, neutral pions promptly decay into two gamma-rays, while the cooling of 
charged pions before their decaying is taken into account. For charged particles, such as 
protons, pions and muons, synchrotron and IC with the full Klein-Nishina cross section are 
included as photon production and cooling processes. The particle decays are simulated with 
the Monte Carlo method adopting the lifetime of pions (muons) as 2.6 x 10 _8 7 7r s (2.2 x 10 _6 7 A1 
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s). The energy fraction of muons at pion decay 7r + (7r - ) — > p^v^ipTv^) is approximated as 
m fj,/ m n ~ 0.76, and the rest of the energy goes to a neutrino. On the other hand, we assume 
that the energy of a muon at its decay fi + (fi~) — > e + u c u^{e~u c u^) will be shared equally by a 
positron (electron), neutrino, and antineutrino. We neglect neutron decay, whose timescale 
is much longer than all the timescales we consider in this paper. 

We also take into account the effect of adiabatic cooling for charged particles with 
the same method in AM11. This effect is controlled by the expansion law of the volume 
V = 4-irR 2 W, where W is the shell width. The volume expansion law affects the escape 
rate for neutral particles. As mentioned in AM11, the escape fraction of photons per unit 
time is c/2W '. We consider several cases for the expansion law in our simulations. 

Assuming the jet opening angle of 9 ]Ct = 10/T, the photon spectrum evolution for 
an observer is calculated in the same manner in AM11. The energy and escape time of 
photons coming from a surface with an angle 9 (measured from the central engine, 8 = 
corresponds to the line of sight) are transformed into energy- and time-bins in observer's 
frame as e = Se'/ (1 + z) (5 is the Doppler factor) and 

t obs = (1 + z) [(1 - & h cos0)rY + Roil - cosd)/c] , (5) 



respectively. When cos# < /3 s h = y 1 — 1/T 2 , photons are coming from the backside of the 
shell with an extra "time delay" due to the shell thickness, 

W'T 

At cx = (1 + z)— [((31 + 1/r 2 ) cos6- p sh ] , (6) 

relative to the emission from the fore side. 

The photon absorption due to extra galactic background light (EBL) is estimated by 
accumulating the opt ical depth between the source and observer with the EBL model of 



Kneiske et al. I ( 120041 ) (best-fit). 



3. Leptonic Models 

3.1. Delayed Inverse Compton 

Inverse Compton emission, which may dominate the GeV energy range, needs low-energy 
seed photons. In the synchrotron self-Compton model (SSC), the seed photons are syn- 
chrotron photons emitted from power-law distributed electrons. The growth timescale of the 



synch rotron photon field may explain the delay timescale of the GeV onset (IBosnjak et al 



20091 ). This possibility was tested by time-dependent simulations in AM11. They showed 
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that a very low magnetic field (e^/e c ~ 10 3 and Rq = 6 x 10 15 cm in their example) is 
required to reproduce the GeV delay. Here, we write the magnetic field 



B' = (7) 



e P cR 2 Y 2 



1/2 / -p \ -1 / d \ -1 / t \ 1/2 



~ V 10 3 / V800; Vl0 15 cm; Vl0 54 er g s-V 

The standard synchrotron model attributes the spectral peak energy £ pea k to the typical 
photon energy emitted from the lowest-energy electrons at injection. The physical mechanism 
that determines the minimum-Lorentz factor 7 C) mm of those electrons is still under discussion. 
In a low magnetic field as given in eq. (jHJ), to emit synchrotron photons in MeV range, a 
very high minimum-Lorentz factor (^> m p /m e ) is required as 



! ' min ~ A/ heBT 



(9) 

where we have omitted the cosmological redshift factor. In this case the typical photon 
energy in the electron rest frame becomes 

^•pcak ^ Tc.min^peak Te,min^peak/r (H) 

- ^UeW VlO" 3 ; V800 ^ VlO^cmy' V 10 54 erg s" 1 J 6H, j 

which largely exceeds m e c 2 , so the Klein-Nishina effect is crucial to emit IC photons. The 
IC photon production is an inefficient process so that the GeV component shows a slow 
growth compared to synchrotron. Soft photons (e <C £ pea k) that can efficiently interact with 
electrons of 7 e , m in are produced by electrons that have already undergone energy losses and 
have 7 e < 7 c ,min- However, as soft photons increases, the newly injected electrons cool via 
mainly IC emission owing to a very low e^/cc- Therefore, the synchrotron component starts 
decreasing while the IC component grows. This enhances the difference in the peak times of 
the MeV and GeV lightcurves. 

For the purposes of comparison with the leptonic calculations in the present paper, we 
summarize here some of the leptonic results in AM11. The lightcurves in Figure 12 in AM11 
show a GeV delay due to the mechanism we mentioned above. However, the difference 
in the peak times is within a pulse timescale, ~ R/cT 2 (the difference in the peak times is 
roughly half of this timescale in the example in AM11), and the GeV and 100 keV lightcurves 
converge in the late stages. The FWHMs of the two lightcurves are almost the same. The 
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example in AM11 adopted an ideal set of parameters for the GeV delay, and it seems difficult 
to produce even larger delays by simple SSC models. We did not find significant GeV delays 
for more conservative SSC parameter sets in our simulations in AM11. 



An external inverse Com pton (EIC) model (IBeloborodov 1120051 ; iToma et al. 1120091 . 12010 



Li 


2010; 


Murase et al. 


2011) 



smaller radius, and they are upscattered by non-thermal electrons in an outer region, where 
some dissipation mechanism, such as internal shocks, inject those electrons. In this model, 
we require two different origins for the MeV and GeV components. A GeV delay larger 
than the typical pulse timescale (R/cT 2 ) appears naturally, although the model requires 
a large number of accelerated electrons, which leads to an energy-budget problem. If the 
numbers of electrons and protons are the same, most of the energy carried by the protons 
remains unreleased. Some GRBs, in which the ma in MeV component is hard to explain 



by the usual synch rotron model, like GRB 090902B (lAbdo et al. 



2009bURvde et al. 112010 



Zhang et al .1120111 ). are interesting applications for this model. In the frame of the outer 
dissipation region, the photon field coming from an inner region is highly beamed along the 
radial direction. The e7-scattering probability becomes highest for head-on collision so that 
the scattered IC-photons tend to propagate with a large angle to the radial axis. As a result, 
the intensity of the scattered photons becomes anisotropic. In this case, the emission from 
off- axis regions (6 > 1/T) contributes considerably. AM 11 showed that the GeV lightcurve 
has a long tail (see Figure 15 in AM11), which is a characteristic feature of this model. 



3.2. Opacity Evolution in Leptonic Model 

As a shell expands, the photon density decreases so that the opacity against 77 absorp- 
tion als o decreases with tim e. This effect may cause a gradual increase of the 77 cut-off 



energy flGranot et al. 1 120081 ) . In this section, we test the evolution of the 77-opacity as a 
possible mechanism of the GeV delay. When the electron injection timescale is comparable 
to or shorter than the shell expansion timescale t' exp = Ro/Tc, the opacity evolution tends to 
cause a negative GeV delay: the GeV emission is terminated earlier than the MeV emission 
as shown in Figure 3 or 8 in AM11. As the electron injection builds up the photon density, 
77-optical depth increases, and GeV photons begins to be absorbed. It was hard to observe 
the effect of the opacity decrease in our simulations with an injection timescale of Ro/Tc. 
Therefore, here we adopt a very long electron injection to see the opacity decay as the shell 
expands. 

We consider a shell expanding from an initial radius Ro = 10 15 cm with a bulk Lorentz 
factor T = 500, and optimize parameters to observe the opacity effects as below. If we take 
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a standard assumption for the shell width as W = Rq/T, the photon escape timescale and 
shell expansion timescale become comparable. In such a case, the opacity evolution may 
be softened by residual photons generated earlier but before escape. So we assume a thin 
shell W = i?o/10r. The electron injection continues until R = 10 16 cm (t' in ^ = 9t' cxp ) with a 
constant rate, and its isotropic-equivalent energy for this pulse is -E e ,pis = 2 x 10 54 erg. The 
injection spectrum is assumed as a cut-off power-law shape, N^i n j(s' e ) oc e'~ Pc exp(— s' e / £ 'eraax) 
for e' c > e' c min , where e' e max is determined by equating the cooling timescale t' cool (synchrotron, 
IC, and SSA are taken into account) and the acceleration timescale 



4c = £^7- (13) 

In this section we assume the electron index p c = 2.5, and "Bohm limit" acceleration as 
£ = 1. In order to avoid contamination of other evolutionary effects, the magnetic field 
B' = 3.9 x 10 4 G, shell width, and minimum injection energy of electrons e' cmin = 2 GeV are 
assumed to be constant during the electron injection. However, the ratio e^/cc, which may 
be approximated as 

e B _ B' 2 f E^pistg^ 1 _ 9R R 2 TB' 2 



e e 8vr \4nR?W T t[ n .J 2£ e , pls 

evolves from 0.04 to 4 during the electron injection. 

Figure Q] shows the spectral evolution in the shell frame. The spectral peak energy 
e' ak ~ 10 keV arising from the parameter set chosen here is somewhat higher than the 
usual values. This is because we are considering a special case similar to GRB 0809 16C, 
which showed a delayed G eV onset with a source-frame peak energy of (1 + 2)e p0 ak,obs ~ 6 



MeV (jAbdo et al. Il2009cl ). In a strict sense, R = Rq in Figure [T] means that a short time 
f /100 has passed after the electron injection started. At this point, the spectrum is not in 
a steady state for photon production, escape, and 77 absorption. During the initial phase 
for R < 1.5Rq, the opacity still increases with the photon density. At R = 1.5Rq the high- 
energy portion of the spectrum can be approximated by a broken power-law with a photon 
index (3 ~ —2.3 below the break energy e' hr and — 3.3 above that , rathe r than an exponential 



cut-off. This break is similar to the prediction by lGranot et al. I ( 120081 ). in which a thin slab 
is considered as a photon source. However, the mechanism for yielding a broken power-law 
is different for our one-zone model with a finite shell width. In a quasi steady state, the 
fraction of photons, whose annihilation timescale t ann is shorter than the escape timescale 
t' esc ~ W'/c, may be written as t'^ nn /t' csc . Then, the spectral shape may be ^o^O^annAeso 
where n' (e') oc e' 13 is the intrinsic spectrum, whose density is determined by the balance 
between continuous photon production and escap e /volume-expansion without the absorption 



effect. The annihilation timescale is t' ann cx ( jLithwick fc Sari Il2001t lAsano fc Takahara 
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Fig. I. — Opacity evolution leptonic model: spectral evolution of photon density in the 
shell frame with the expanding radius R (see the text in §3.2p . Initially the photon density 
increases as electrons are injected (R < 1.5Rq), then the shell expansion and photon escape 
make the density begin to decrease. In this stage the break energy e' hr ~ 10 7 eV due to 
77-absorption gradually grows as the shell expands. 



20031 ). while t' csc is independent of the photon energy. Therefore, n'{e') oc e' 2 P +1 above the 
spectral break at e' hT ~ 10 7 eV. The value (3 ~ —2.3 implies that the index should be —3.6. 
We should note that photons above s' > 10 7 eV interact with photons below e' pcak . In our 
example in Figure HJ the index of the photon distribution for one order of magnitude just 
below e' peak is a ~ —1.3. In such cases, the spectrum should be n'(s') cx £ , l 3+a+1 ~ e'~ 3 - . 
Our result seems consistent with those estimates. 

In the later stage (R > 2Rq), the 77-attenuation gradually becomes inefficient, and this 
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causes the photon spectral break energy to increase with time. In this stage the photon 
distribution may be in the quasi steady state. Therefore, the spectra for R = R (photon- 
density growing stage) and R = 5Rq are different in spite of their similar photon densities. 
The deviation from a s imple power-law seen around 1 eV is due to electron heating via SSA 
fiGhisellini et al. Ill988i AM11). 



ef(e) [erg/cm /s] 




£[eV] 



Fig. 2. — Opacity evolution leptonic model: flux evolution for an observer with z 
(see the text in §3.2j) . 



4.35 



Assuming a source redshift z = 4.35 such as that of as GRB 080916C, we show the evolu 
tion of the spectrum for an observer in Figure El Given a luminosity Lh above 

£pen,k ^ m a case 

of the usual short electron- injec tion (t[ n ^ < Rq/Tc), an analytical estimates (ILithwick fc Sari 
200ll ; lAsano fc Takahara 1120031 ) gives us the 77 cut-off /break energy as 



'77 



2 4 

— ~ — r P+ 1 

£pcak 



lQTTE peak C 2 8t 



oc St w+i} 



(15) 
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where the dimensionless function F(/3) ~ 10~ 2 , and St = R/T 2 c is the variability timescale. 
For our continuous injection case, that may be interpreted as e 77 oc t^ s with /3 = —2.3. By 
comparing spectra for t G b s = 0.31 s and 3.1 s, we can see that the spectral behavior is close 
to this analytical approximation; the break energy grows by a factor of six or more. 

Normalized Flux 




Fig. 3. — Opacity evolution leptonic model lightcurve (see the text in §3.2p . The source 
redshift is assumed to be 4.35. The lines labeled with "GeV" or "MeV" are plotted based 
on the spectral flux, while the line labeled with ">GeV" is a curve based on the photon 
number integrated above GeV. 

However, the spectral shape does not show a sharp cut-off or break. The deviation 
from the extrapolation of the low-energy spectral shape is gradual. As a result, the spectral 
flux at GeV (see "GeV" in Figure |3]) shows a similar evolution to the flux at MeV. The 
MeV lightcurve has a plateau reflecting the constant electron injection rate. The peak time 



13 



of the GeV lightcurve delays relative to the MeV arising time, but the GeV flux changes 
slowly during the MeV-plateau phase. However, for low photon statistics, the lightcurve 
would be obtained by accumulating photon counts above GeV. As shown in Figure [31 the 
lightcurve based on the accumulated photon counts shows a distinct delay relative to the 
MeV lightcurve. 

We have considered an artificial setup to see a GeV delay due to the opacity evolution. 
Even in this extreme case, the spectral evolution is gradual, and it seems different from the 
observed "sudden" onset of GeV emission. 



4. Hadronic Models 



The electromagnetic cascade triggered by photop ion production can make GeV extra 



compo ne nts a s shown in our series of GRB studies (lAsano fc Inoue I l2007t lAsano et al. 



2009al lbl. [2010l l The electron cooling timescale is much shorter than the timescales of pro- 
ton acceleration and photopion production. Those timescales and continuous photomeson 
production after the end of electron injection may cause a delayed onset of the GeV emission. 

In this section, we inject electrons and protons with the same timescale t[ n ^ = R /Yc. 
The injection spectra have the same shape as that for electrons in §3.21 The proton injection 
i ndex is assumed to be p p = 2.0, which is appropriate to make GRBs contribute to UHECRs 
(IWaxman fc Bahcall 1119981 1. while that for electrons is taken to be p ? = 3 . to reproduce 
the ty pical photon index (3 ~ —2.5 measured by BATSE (IPreece et al. Il2000t iKaneko et al. 
20061 ). Although the injection indices for protons and electro ns are different, the en ergy scales 
we consider are largely different for electrons and protons (see lAsano et al. Il2009bl . for further 
discussion). We adopt the parameter £ = 1 for both electrons and protons (see eq. f JT3|) ). 
The maximum energy of protons is determined by equating the cooling and acceleration 
timescales or by the condition that the Larmor radius should be shorter than the shell 
width. We take into account proton synchrotron and photopion production to estimate the 
cooling timescale. While the proton maximum energy at injection is controlled by the time 
step for injection, we phenomenologically carry out succeeding proton acceleration following 
the acceleration timescale and index. The minimum energy of the protons is fixed at 3 GeV. 

We find that the shell width expansion enhances the GeV delay in our simulations. 
Below, we show examples with an expanding width as W = R/Y. During the particle 
injection, we neglect adiabatic cooling and decay of the magnetic field. After the particle 
injection, the magnetic field is assumed to decay as B' oc R~ 2 . 



In this paper the MeV photon peak is assumed to be produced by accelerated electrons 
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i njected in th e emission region, which differs from the previous assumption in lAsano et al. 



( I2009all2010l ). where an ad-hoc Band- type photon component matching the observed spectra 
was postulated. Therefore, differences from the preceding simulations arise not only because 
here we consider the time-dependence, but also because some of the basic input assumptions 
are different. One issue here is that the synchrotron emission from the primarily accelerated 
electrons yields a soft photon index a ~ —1.5, which is difficult for reproducing the observed 
spectral indices. This is because the effective energy-distribution for cooled electrons below 
7e,min becomes n e {^ e ) oc 7~ 2 , when the synchrotron emission is the dominant cooling process 
(7 e oc 7g). In order to resolve this contradiction, another mechanism for the Band component 
i s required, such as a pho t osphere model or co ntinuous particle acceleration by turbulence 



f lAsano fc Terasawa 1 12009c iMurase et al. 1 120121 ) . However, since our emphasis here is more 



on the GeV components, for simplicity we consider here the standard synchrotron model. 



4.1. Hadronic "Moderate" Case 

Here we adopt a typical GRB luminosity for one pulse in a GRB, as opposed to the 
high values deduced for LAT bursts. The model parameters are V = 800, Rq = 10 14 cm, 
i?e,pis — 5.0 x 10 50 erg, and an isotropic-equivalent energy of protons -E PiP i s = 6.1 x 10 51 erg 
(e p /e c ~ 12)0 including the acceleration effect after injection. The initial magnetic field B' 
is set through the parameter ratio at R = 2R Q as 

6r ARnB!? 

— ~ = 3.0, (16) 

-£-'e,pls 

which results in B' = 1.9 x 10 4 G. The peak energy e pca k is adjusted to be ~ 100 keV by 
adopting e' p min = 380 MeV. 

In Figures H] and O the main particle spectra are plotted in the shell frame, omitting for 
clarity the spectra of the electron/positron, pion, muon, and neutrino components. In the 
early phase (R < 2R ), the proton density is increasing by proton injection, and their maxi- 
mum energy is growing with the radius R obeying the proton acceleration timescale. Then, 
protons undergo adiabatic cooling, while neutrons keep their energies until they escape or 
interact with photons. The volume expansion (V oc R 3 ) reduces the particle number den- 
sities, and particle escape also affects the densities of neutrons and photons. The evolution 
of the neutron-to-proton ratio from 2i? to 3Rq indicates that pion (and neutron as well) 
production continues even after the end of particle injection. Therefore, the electromagnetic 
cascade due to photopion process continues as well for t' > t- nj . 



J An alternate notation in some papers is l// e ~ ep/ee 
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Fig. 4. — Hadronic model, "moderate" case: spectral evolution of the photons, protons, and 
neutrons in the shell frame with the expanding radius R (see the text in §4.11) . As particle 
injection proceeds, the photon and proton densities build up. Then, for R > 1.5Rq, the 
density starts to decrease due to the effects of the shell expansion and photon escape. The 
electromagnetic cascade triggered by photopion production gradually enhances the height of 
the second peak (e' ~ 10 8 eV) relative to that of the main peak (e' ~ 10 2 eV) in the photon 
spectrum. 

In the photon spectra, we can see a 77 cut-off at ~ 10 8 eV and an SSA (synchrotron 
self-absorption) signature at ~ 0.1 eV. The second peak of the gamma-ray spectra at ~ 10 8 
eV is due both to SSC by the primarily accelerated electrons and to synchrotron emission 
by secondary electrons/positrons. The cascade due to photopion process raises the second 
peak relative to the synchrotron-first peak even for R > 2R . At R = 10i?0) the highest 
energy region (> 10 16 eV) is dominated by neutrons, because protons have been adiabatically 
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Fig. 5. — Hadronic model, "moderate" case: same as Figure H] but for later stage. Neutrons, 
which, being neutral, are not subject to adiabatic cooling, dominate the high energy region 
above 10 15 eV. 

cooled. 

Assuming z = 1.0, we plot the spectral flux evolution of gamma-rays for an observer in 
Figure [6j It is clearly seen that the spectral shape above £ pca k becomes harder with time, 
and the fraction of the GeV-10 GeV component grows. The softening in the lower-energy 
region (e < e pca k) is also caused by the hadronic cascade. The time-integrated spectrum is 
shown in Figured in which we compare this hadronic model with the same (leptonic) model 
without proton injection. Both models show a second peak at 10 10 eV, in the leptonic case 
this being due to the SSC emission, while in the hadronic it is due to the cascade. Thus, it is 
hard to distinguish these model just from their spectral shape. However, as shown in Figure 
[HI the lightcurves for the hadronic model (solid lines) may be distinguishable from those for 
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8 [eV] 

Fig. 6. — Hadronic model, "moderate" case: flux evolution (see the text in §4. 1[) for an 
observer with z = 1. 

the leptonic models (dashed lines). One sees that the 100 keV and 10 GeV lightcurves for 
the leptonic model (dashed lines) have almost the same shape, and the GeV delay does not 
appear in this leptonic case. In order to produce a GeV delay in these leptonic models, we 
may need a smaller magnetic field, such as e^/ec < 10 -3 or B' < 100 G, as shown in AM11. 
On the other hand, for the "moderate" hadronic model, the 10 GeV or 30 GeV lightcurves 
show a distinct (even if short) delay relative to the 100 keV lightcurve. Even so, the peak 
time for the 10 (30) GeV-lightcurve is earlier than the peak time for the neutrino lightcurve. 
This may be because the SSC contribution from primary electrons is not negligible in the 
10 (30) GeV-lightcurve. 

In this example, the isotropic-equivalent energies of escaped neutrinos and neutrons are 
4.4 x 10 49 erg, and 2.2 x 10 50 erg, respectively, per pulse. For reference, let us roughly 
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Fig. 7. — Fluences for the hadronic model and purely leptonic model (e p = 0) with the 
"moderate" parameters (see the text in §4.11) for an observer with z = 1. The dashed line is 
the spectrum without the absorption effect of EBL. 



estimate diffuse neutrino background and the contribution of the escaped neutrons to the 
UHECR flux based on this specific parameter set. Here we consider multiple pulses with 
the identical parameter set for a GRB. The pulse number N = 30 with the "moderate" 
parameter set yields the total isotropic-equivalent gamma-ray energy of Ei SO = 2.0 x 10 52 
erg. This energy scale can be used as a typical example. The local GRB rate ~ 1 Gpc~ 3 yr" 1 
with the average energy E iso = 2.0 x 10 52 erg corresponds to the gamma-ray energy release 
rate of 2 x 10 43 erg Mpc ~ 3 yr -1 . The redshift z = 1 in this example is lower than the 



observed average redshift (IJakobsson et al. 1 120061 ) . so the total fluence 7.4 x 10 6 erg cm 



per burst is also slightly larger than the typical GRBs. Si mply accumulating this identical 
GRB at z = 1 with N = 30 and the GRB rate of 667 yr" 1 jAbbasi et al. II2012I ). we plot the 
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Normalized Flux 




Fig. 8. — Lightcurves for the "moderate" case with z = 1. See the text in §4.11 The 
designation of the labels is the same as that in Figure [3J The solid lines correspond to 
the hadronic model, while the dashed lines are for the same leptonic model without proton 
injection. The neutrino flux evolution at 10 16 eV is also plotted. 



spectra of the neutrino background in Figure Here we have not distinguished between the 
neutrino species, so the spectrum consists of neutrinos originating from both pion-decay and 
muon-decay. Even though the adopted gamma-ray ffuence is larger than the a verage one, 
the neutrino fluence is lower than the IceCube limit for (lAbbasi et al. II2012I ). This may 
be due to the large Rq and V in this model. Furthermore, note that the dominant energy 
region for this mo del is 10 17 -10 1 ^ eV, w hich is well above the energy region constrained with 
IceCube. Both in lAbbasi et al. I (120121 ) and here, the neutrino spectrum is consistent with 
the commonly considered case where neutrinos from pion-decay dominate, in which the lower 



-20 - 




Fig. 9. — Spectra of the diffuse neutrons and neutrinos (without the GZK effect) adopting a 
model corresponding to the "moderate" case (e p /e c ~12), from z = 1 (full lines, accumulating 
667 GRBs yr _1 ), and from lOOMpc (dashed line, for neutrons with a local rate 2 GRBs 



Gpc 3 yr 1 ); see the text in §4.11 The upper limit for der ived by IceCube ( Abbasi et al 



20121 ) and the UHECR flux determined by the Auger team ([Abraham et al. 
plotted. 



201011 are also 



energy break of the neutrino spectrum (IWaxman fc Bahcall 1119971 ) is given by 



8.6 x 10 16 



-1 



V lOOkeV/ V800 / 6 



(17) 



where we take 5 = 2T as the Doppler factor, and K pi = 0.2. As shown in lAsano et al. 
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f l2009a| ). a large Lorentz factor or a large radius Rq is required in order to have GeV photons 
which escape unabsorbed from the source, and this leads to a lower photopion efficiency 
and a harder neutrino spectrum. Therefore, such models which result in GeV emission 
can avoid the GRB- neutrino upper limit constraints given by IceCube. Also, the report in 



Abbasi et al. I (120121 ) excludes only GRB models with T < 400. 



With the same method as for the neutrinos, we also plot the spectrum of the escaped 
neutrons without the GZK effect in Figure O The neutron energy range extends as far as 
the UHECR energy region, and the contribution from those 6 67 GRBs in this simple model 



seems close to the UHECR spectrum obtained with Auger ([Abraham et al. 1 120 10h . The 



isotropic-equivalent energy of the neutrons escaped per burst (assuming N = 30 pulses) is 
6.6 x 10 51 erg, comparable to the energy emitted per burst (30 pulses) in 7-rays, 2 x 10 52 erg. 
However, at 10 20 eV the GZK effect is crucial so that we should count only GRBs within 
the GZK horizon (~ 100 Mpc). The GRB rate within 100 Mpc sphere (~ the GZK horizon) 
becom es 8.4 x 10~ 3 GRBs yr _1 for the optimistic rate of 2 Gpc -3 yr -1 ( Wanderman &: Piran 



20 lOn . Here we simply assume that GRBs occur with this rate at a distance of 100 Mpc, 



and we plot the UHECR flux from those GRBs as the dashed line in Figure O Note that 
even though we adopt e p /e c — 12. the energy fraction of escaped neutrons to gamma-rays 
is ~ 1/3, due to the larger radius and T leading to a modest pj rate. The UHECR flux 
level at 10 20 eV is lower than the observed flux by a factor of ~ 40, which implies that 
ten times the gamma-ray energy is required to be emitted as cosmic rays to make GRBs 
dominant sources of U HECRs. This is reco nfirmation of the difficulty in GRB-UHECR 
scenario pointed out in lEichler et al. I (120101 ). Here we have taken the most pessimistic 
assumption for UHECR production, namely, that the high-energy protons cannot escape 
until they have cooled adiabatically. As a result, the escaped-neutron energy is a few percent 
of the injected energy of protons. Some mechanism to make protons escape may be required 
to produce UHECRs from GRBs. Of course, we also have to include the lum inosity function, 



redshift distribution, and the uncertainties discussed in IWaxman I (120101 ). Such detailed 



estimate of neutrinos and UHECR with our numerical simulations will be addressed in our 
future work. 



4.2. Hadronic "Fermi-LAT" Case 

Here we adopt a more specific parameter set that mimics a pulse of the bright GRBs 
observed with Fermi-LAT . The model parameters are T = 600, Rq = 1.3 x 10 16 cm, -Ec.pis = 
2.0 x 10 54 erg, and i? P;P i s = 4.8 x 10 55 erg (e p /e e — 24) including the acceleration effect after 
injection. The initial magnetic field B' = 830 G is determined by the same method in §4.11 
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as es/e c = 3.0 at R = 2R . The peak energy £ pe ak is adjusted to be ~ 1 MeV by adopting 
<mm = 13 GeV. 



e' 2 n'(e') [erg/cm 3 ] 




Fig. 10. — Hadronic "Fermi-LAT" case: spectral evolutions for photons, protons and neu- 
trons in the shell frame with the expanding radius R (see the text in §4.2p . The bump in 
the photon spectra at 10 7 -10 10 eV is due purely to hadronic cascades. 

The evolution of the particle spectra in the shell frame (Figure [TOl) is similar to that 
for the "moderate" case in Figure HI However, the gamma-ray spectral bump above 10 7 
eV is here attributed purely to hadrons. Although the ratio e#/e e is the same as that 
for the "moderate" case, the Klein-Nishina effect is crucial because of the very high e' e min . 
Therefore, the IC emission does not yield a spectral bump in this case. As the photon density 
increases, photons due to electromagnetic cascade triggered by photopion production appear 
above 10 7 eV. 
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Assuming z = 4.35, we plot the observer-frame spectral evolution in Figure [TT1 Initially 
(t < 1.6 s) the spectral shape is close to the simple Band function, and no spectral bump 
appears in the GeV energy range. In the later phase, the hadronic component appears, and 
the entire shape gradually evolves into flat on e. This evolution is similar to the spectral 
behavior in GRB 080916C (jAbdo et al. 112009a ). The effect of the hadronic cascade on the 
spectral flatness is well shown in the time- integrated spectra in Figure [121 In contrast to this 
two-component hadronic behavior, the purely leptonic case for the same parameters leads to 
a pure Band-type time-integrated spectrum. The spectral shape is almost unchanged during 
the emission in the leptonic model. 



The lightcurve in Figure [13] shows a significant delay of the 10 GeV-lig ht curve. The 
delay timescale (~ 6 s) is comparable to the delay seen in GRB 080916C (jAbdo et al. 
2009a ). In this case SSC is negligible, so the GeV spectral-bump is purely due to hadronic 
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Fig. 12. — Fluences for the hadronic model and purely leptonic model (e p = 0) with the 
"Fermi-LAT" parameters (see the text in §4.21) for an observer with z = 4.35. The dashed 
line is the spectrum without the absorption effect of the EBL. 

cascade. Unlike in the "moderate case" , the peak time of the 10 GeV photon lightcurve is here 
relatively close to the peak time for the neutrinos. Thus, the lack of SSC effects (moderated 
by the Klein-Nishina effect here) is of key importance in emphasizing the delayed onset in 
hadronic models. The 10 GeV- light curve shows a significantly broader FWHM (~ 14 s) than 
the FWHM of the MeV-lightcurve (~ 11 s). In the same model without proton injection 
(dashed lines), the high-energy lightcurve is rather narrower than the MeV lightcurve. The 
onsets of the GeV and MeV lightcurves are almost simultaneous. Thus, for this parameter 
set, the difference between hadronic and leptonic models is prominent. 

Here we have adopted a very large energy per pulse. In such cases the number of pulses 
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Fig. 13. — Lightcurves for the "Fermi-LAT" case with z = 4.35. See the text in §4.21 The 
designation of the labels is the same as that in Figure EJ full lines are for the hadronic model, 
dashed are for the purely leptonic one. In the hadronic model, since the 10-30 GeV photons 
originate mainly from the hadronic cascades, the peak times for those lightcurves and for 
the neutrino curve almost coincide with each other. 

may be a few, and the expected rate of such bright Fermi-LAT GRBs (in this example the 
ffuence is 7.5 x 10~ 5 erg cm -2 ) may be less than 10 per year. Thus, in the same manner 
as in Figure [9] adopting N = 1 and the rate 10 yr~ x , we plot the spectra of the neutrino 
and neutron due to bright Fermi-LAT GRBs in Figure HH The spectrum of the escaped 
neutrons ranges as far as the UHECR energy scale in this model also. The contribution to 
the diffuse neutrino background and to the integrated flux of escaping neutrons are smaller 
than in the "moderate case", owing to the low GRB rate with such parameters. However, 
the spectral flux of the neutrons at 10 20 eV is comparable to that in the "moderate case", 
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Fig. 14. — Spectra of the diffuse neutrons and neutrinos adopting a model with the "Fermi- 
LAT" hadronic cascade case (accumulating 10 GRBs yr -1 ). See the text in §4.21 

although such neutrons lose their energy by the GZK process. (Note that here we assumed 
z = 4.35, and the equivalent rate from within 100 Mpc would be much smaller than that 
for the "moderate" hadronic case parameters). The isotropic-equivalent energy of neutrons 
per pulse is 8.7 x 10 53 erg, which is only 1.8 % of the injected energy of protons, while the 
emitted isotropic-equivalent gamma-ray energy per pulse is 2.8 x 10 54 erg. 

The isotropic-equivalent energy of the neutrinos per pulse is 1.8 x 10 53 erg. The break 
energy at ~ 10 17 eV in the neutrino spectrum is in apparent contradiction with equation 
(TIB"]) for this parameter set. The reason is that in this case the pion production efficiency 
is not high enough, even for protons of 7 p = 7 P)th , the efficiency increasing as oc 7p 5 , as 
mentioned in $2] The neutron spectrum in the "Fermi LAT" case in Fig. [14] has a sharper 
peak at ~ 3 x 10 19 eV than the "moderate" model of Fig. O This indicates that only very 
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high-energy protons near the cut-off energy produce pions efficiently. The high-energy bump 
at ~ 3 x 10 18 eV in the neutrino spectrum is attributed to pion-decays from such high-energy 
protons/neutrons, while the low-energy break is attributed to the muon-decay component. 
Thus, also for this parameter set, the neutrino energy range is well above the energy range 
constrained by IceCube. The neutrino fluence (~ 1 0~ 6 erg cm" 2 ) from th is pulse corresponds 



to ~ 1/300 of the IceCube limit (~ 0.2 GeV cm- 2 . lAbbasi et al. 112012 ). 



4.3. Proton Synchrotron 



A different hadron ic model is the proton synchrotron model (jVietri 111995c ; iTotani 111998 



Razzaque et al. 1 120101 ) . To test this case with the same Lorentz factor and initial radius as 
those in §4.21 we need to use extreme values for some of the other parameters, in order to 
ensure that proton synchrotron contributes to GeV energy range. Thus, we adopt Ee.pis = 
4.0 x 10 54 erg, and E PtPis = 2.0 x 10 56 erg (e p /e c = 50), and the initial magnetic field B' Q = 4800 
G (es/e c = 50 at R = 2R ). The minimum energy of electrons is set as e' cmin = 5.3 GeV. 
If the acceleration parameter £ for protons is unity, the maximum energy of protons in this 
strong magnetic field becomes too large and cascade emission from photopion production 
overwhelms the proton synchrotron emission. This is because the soft photon spectrum 
a ~ —1.5 implies / P7 oc 7p' 5 . In order to suppress photopion production, we degrade the 
maximum energy of protons by adopting £ = 100. 

The resultant spectral evolution for an observer at z = 4.35 is shown in Figure [TBI Here 
the second peak of the photon spectra is due to proton synchrotron radiation. It is seen that 
proton synchrotron can produce harder spectra than the cascade emission seen in Figure [TT1 
Even in this case, however, we cannot neglect photopion production. In the later phases, the 
spectrum becomes flat owing to the cascade triggered by photopions. The delay time of the 
10 GeV-lightcurve (Figure ITB]) is similar to that for the "Fermi-LAT" case. This is because 
the parameters T and R are common and the GeV bump is due to purely hadronic effects 
in both models. However, the FWHM of the 10 GeV-lightcurve is narrower here and the 
decay timescale of the 10 GeV emission seems almost the same as the timescale in the MeV 
emission, while the 30 GeV-lightcurve is narrower than the MeV lightcurve. This depends on 
the redshift, or on model parameters (r etc.). However, for the proton synchrotron model, 
the pulse decay timescale roughly corresponds to the synchrotron cooling timescale, so that 
the decrease of the pulse timescale with photon energy may be a characteristic feature of 
such models. Although the degeneracy in the parameters is not easy to resolve, having an 
energy threshold as low as possible could help in distinguishing the previous "Fermi-LAT" 
hadronic cascade and this "Fermi-LAT" proton synchrotron model through their lightcurve 
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Fig. 15. — Proton synchrotron model: flux evolution (see the text in §4.3}) for an observer 
with z = 4.35. 

differences. 

The emitted energy as gamma-rays, neutrons, and neutrinos per pulse are 9.0 x 10 54 
erg, 3.2 x 10 54 erg, and 8.0 x 10 53 erg, respectively. In Figure [T71 we plot the spectra of the 
escaped neutrons and neutrinos in the same manner (N = 1, 10 GRBs yr^ 1 ) as in Figure 
IT41 The maximum energy of neutrons is suppressed by the large £ (see equation (I13p ). The 
neutron and neutrino spectral shape are similar to those in Figure [141 respectively, although 
the typical energy ranges become lower. This is because the critical parameters for photopion 
production (R , T, and Ee.pis) are almost the same in the two cases. 
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Fig. 16. — Proton synchrotron model: lightcurves observed from z = 4.35. See the text in 
§4.31 The designation of the labels is the same as that in Figure |8j 

5. Summary and Discussion 

We have simulated the time evolution of the particle energy distributions in a relativis- 
tically expanding shell with parameters corresponding to that of GRBs observed by Fermi, 
and discussed the spectral evolution of the radiation seen in the observer frame. If we adopt 
a very weak magnetic field in the shell in a one- zone leptonic model, the slow evolution of the 
synchrotron self-Compton (SSC) emission due to the Klein-Nishina effect results in a delayed 
onset of the GeV emission, by amounts comparable to those observed by Fermi. However, 
in such leptonic models the FWHM of the GeV lightcurve is almost the same as that of 
the 0.1-lMeV lightcurves, although the present-day low photon statistics in the GeV range 
makes it hard to compare the FWHMs. On the other hand, an external inverse Compton 
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Fig. 17. — Spectra of the diffuse neutrons and neutrinos adopting a model with the "Fermi- 
LAT" proton synchrotron case. See the text in §4.31 

(EIC) model, which is essentially a leptonic two-zone model, can reproduce the GeV-delay 
naturally. A long tail for the GeV lightcurve, as shown in AM11, is possible evidence for this 
model. T his can partially contr ibute to the extended GeV emission seen in several GRBs 



(see, e.g., iGhisellini et al. 1 120101 ) . We have also tested the 77-opacity evolution effect with 
our code. At least in our one-zone formulation, an electron injection timescale longer than 
the expansion timescale is required in order to observe the opacity damping. The resultant 
lightcurves have a characteristic shape, and the growth of the cutoff energy is more gradual 
than in the current sample of Fermi-LAT bursts, which so far do not show a similar spectral 
evolution. Thus, 77-opacity evolution cannot be the reason for the delay. 



The hadronic models studied here include both cascade models and a proton synchrotron 
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model, both of which are able to reproduce the delayed onset with a delay timescale close 
to the pulse timescale (~ Rq/cT 2 ). The delay is due to the long acceleration timescale of 
protons and continuous photopion production after the end of the particle injection. The 
wider FWHM for the GeV lightcurve than for then MeV lightcurve could be a signature of 
the hadronic cascade. If the Klein-Nishina effect prevents IC emissions, the delay due to 
hadronic cascade becomes more dominant. The amounts of escaped UHE neutrons in our 
examples are comparable to the gamma-ray energy. Since we have adopted large T and R 
values in order to simulate GRBs where GeV photons are able escape from the source, the 
resultant neutrino spectra are hard enough to av oid the current flux limi t constraints from 



IceCube. In addition, as the Fermi team indicates ( lAckermann et al.ll2012[ ). GRBs with extra 
components in the GeV band may be a small fraction of the GRB population. Therefore, 
even if GRBs accelerate UHE particles, the typical GRB parameters are in a range which 
implies a sufficiently low neutrino and GeV gamma-ray production efficiency Our hadronic 
model also predicts a delayed onset of the neutrino emission, which is more pronounced 
than the corresponding GeV photon delay. If neutrinos are eventually observed, this may 
be another point which could be tested by future neutrino observations. The energy budget 
required for being UHECR sources appears insufficient in our examples. In future work, we 
plan to include the redshift evolution, luminosity function, and other details to further test 
the viability of the GRB-UHECR scenario. 

The open problem of the low-energy photon spectral index a is not addressed in this 
paper. Here we inject the electrons in the usual manner; the minimum energy of the electrons 
at injection e' e min being a free parameter^. The synchrotron cooling results in a soft effective 
electron distribution (n' c (e' c ) oc e'~ 2 ), so a < —1.5 is unavoidable in this type of simulations. 
Spectral slopes such as these, while softer than the average a ~ — 1, are present in the 
BATSE data base. Such softer slopes have the result of enhancing the photopion production 
efficiency, as well as the resultant flux of escaped neutrons. Note that here the energy 
fraction in protons required to generate a GeV extra component is not so high, e p /e c = 10- 
25 (~ f~ x , in alternate notation), compared with the case of GRB 090510 where e P /e e > 100 



was required to reconcile with the observed index a = —0.58 ( lAsano et al. Il2009af ). Another 
consequence is that the soft spectra make it difficult to find a low-energy spectral excess 
in the X-ray region. In our hadronic examples, the low-energy excess appears in the later 
stages, when the power-law component due to the hadronic cascade becomes prominent or 
dominant over the entire spectrum. On the other hand, in order to have a self-consistent 
one-zone synchrotron model whether leptonic or hadronic, some mechanism to harden the 



2 As is well known, e' e min is written by the phenomenological parameters: e and the number fraction of 
the accelerated electrons, both of which are not well constrained so far. 
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spectrum (e.g., lAsano fc Terasawa 1 120091 ; iMurase et al. Il2012l ) would need to be included, 
which could change the lightcurves or the UHECR production efficiency. 

In our simulations the minimum-Lorentz factor of electrons at injection 7^ min is assumed, 
for simplicity, to be constant during the injection. Thus, our simulations do not show a 
significant £ pea k-evolution, as seen in Figure [HI whereas an increase of £ pca k around the 
time of the GeV onset has been reported in GRB 080916C and 090902B. To reproduce such 
a £ pcak -evolution in our model, a growth of j' p miT1 during the injection would need to be 
introduced. Ho wever, in GRB 090926A the £ pea k does not change drastically around the 
GeV onset time (lAckermann et al.ll201ll ; IZhang et al.ll201ll ). Given this variety of spectral 
evolution behaviors, at this stage one is unfortunately left without a definite guideline for 
including a temporal evolution of the model parameters. 

In summary the GeV delays can be explained either by a one-zone leptonic model with 
very low magnetic field and high luminosity (AM11), or by normal parameter range one- zone 
hadronic models (this wo rk), or also by a n o rmal magnetic field and luminosity two- zone 



leptonic model via EIC (IToma et al. 1 12009k |2010| ). However, the lightcurve behavior of 



these models differ substantially from each other. The qualitative tendency of the multi- 
GeV lightcurves discussed here can be tested with future atmospheric Cherenkov telescopes 
such as CTA. For this purpose, it would be desirable that the low-energy threshold of those 
telescopes should be as low as possible, in order to ensure good photon statistics as well 
as the capability to observe the spectral evolution below the 77 cutoff energy due to EBL. 
Such an ideal instrument would provide the capabilities for significant breakthroughs in GRB 
physics. 
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